clear;
T=6;
N=100000;
g=9.8;
l=0.1;
dt=T/N;%diataT
t=0:dt:T;%时间
y0=2.5*pi;%初始条件
y1=zeros(1,length(t));%建立一个全是0的向量
y1(1)=y0;%将初始条件放入
v0=0.5;%初始条件
y2=zeros(1,length(t));%建立一个全是0的向量
y2(1)=v0;%将初始条件放入
for n=1:length(t)-1
    y1(n+1)=y1(n)+y2(n)*dt;
    y2(n+1)=y2(n)+(-g/l*sin(y1(n)))*dt;
end
%plot(t,y1,t,y2);
plot(y1,y2)
axis equal
%利用holdon更改划分，判断怎么样的精度合适
hold on